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CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION OF AEROELASTIC EFFECTS 
ON THE LATERAL CONTROL OF SWEPT AND UNSWEPT WINGS}! 


SUMMARY 


Charts and approximate formulas are presented for the estima- 
tion of static aeroelastic effects on the spanwise lift distribution, 
rolling-moment coefficient, and rate of roll due to the deflection 
of ailerons on swept and unswept wings at subsonic and super- 
sonic speeds. Some design considerations brought out by the 
results of this report are discussed. 

This report treats the lateral-control case in а manner similar 
to that employed in NACA Report 1140 for the symmetric- 
flight case and is intended to be used in conjunction with NACA 
Report 1140 and the charts and formulas presented therein. 


INTRODUCTION 


The lateral control and maneuverability of a wing are 
Important design considerations. These characteristics may 
be affected to а significant extent by aeroelastic action, 
particularly at high dynamic pressures and in the case of 
thin wings, swept wings, aud wings designed for low wing 
loadings, because the operation of ailerons and spoilers 
usually creates aerodynamic forces which deform the wing. 

As a result of these deformations, the angles of attack 
along the span often change in such a manner as to produce 
lifts which oppose the rolling moment of the aileron or 
spoiler; furthermore, these lifts cause additional deforma- 
tions which may again reduce the rolling moment, and so on, 
until equilibrium is reached. Wing flexibility may thus 
cause & serious loss in the control power; in fact, if the 
dynamic pressure of the airstream is sufficiently high, the 
aileron rolling moment may be completely nullified. The 
speed and dynamic pressure at this condition are often re- 
ferred to as the aileron reversal speed and reversal dynamic 
pressure, because at higher dynamic pressures the controls 
would have to be reversed in order to roll the airplane. 
When wing flexibility causes a loss in lateral control, there is 
also usually a loss in the rolling maneuverability, which may 
be expressed as the wing-tip helix angle due to rolling and is 
affected by changes in both the control power and the damp- 
ing in roll. These aeroelastic effects on the lateral control 
and maneuverability have to be чып into account in a 
design of & wing. 

Several methods аге available for calculating diens 
effects (ref. 1, for instance), but, since these effects depend 
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on the structural characteristics of the wing, which are not 
accurately known in advance of its design, the relatively 
large amount of time required for even the most efficient 
of these methods militates against their use in connection 
with preliminary design calculations. А need exists, there- 
fore, for means of estimating quickly some of the more 
important aeroelastic effects on lateral control with an 
accuracy that is sufficient for preliminary design purposes. 

The related problem of estimating static aeroelastic 
effects on the magnitude and spanwise distribution of the 
lift in symmetric flight has been treated by the charts and 
approximate formulas presented in reference 2. The pres- 
ent report consists of an extension of the analysis of refer- 
ence 2 to the lateral-control case. Inasmuch as the static 
aeroelastic equations are linear, the results presented in 
the two reports may be superimposed. Included in the 
present report are approximate formulas for the estimation 
of the static aeroelastic effects on the spanwise lift distri- 
bution, rolling moment, and rate of roll due to aileron deflec- 
tions on swept and unswept wings at subsonic and supersonic ` 
speeds. Also presented are summary charts which 
indicate whether a given design is likely to be affected by 
losses in lateral control. By means of these charts and 
approximate formulas as well as those of reference 2, the 
conventional procedure of designing a wing on the basis of 
certain strength criteria, checking ic for aeroelastic phe- 
nomena, and then reinforcing 16, when necessary, to meet 
the stiffness requirements imposed by these phenomena can 
often be simplified greatly, inasmuch as the effect of some of 
these phenomena,can be estimated in advance of design. 

In order to keep the length of the report to & minimum 
and to avoid a repetition of much of the material presented 
in reference 2, the present report has been written in such & 
manner as to facilitate its joint use with reference 2 rather 
than to make it entirely self-contained. The use of the charts 
and approximate formulas presented herein is described and 
the limitations of the charts and the light they shed on some 
design problems are discussed. A numerical example is 
included to illustrate the use of the approximate formulas 
of this report. А brief description of the calculations on 
which the charts апа approximate formulas are based is 
contained in the appendix to supplement the more detailed 
derivations in references 1 and 2. 


! Previously released as NAOA TN 2747, “Charts and Approximate Formulas for the Estimation of Aercelastic Effects on the Lateral соно of Swept and Unswept Wings” by 


Kenneth A. Foss and Franklin W. Diederich, 1952. 
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SYMBOLS 

aspect ratio, 02/5 

swept-span aspect ratio, .А/соз% А 

location of section aerodynamic center 
measured from leading edge, fraction 
of chord 

wing span, in 

chord (measured perpendicular ia elastic 
axis), in 

aileron chord, in. 

section lift-curve slope per radian 

effective wing lift-curve slope per radian 

rolling-moment coefficient 

damping-in-roll derivative | 

rolling-moment coefficient due to aileron 
deflection 

location of chordwise center of pressure 
of lift produced by aileron deflection 
measured from leading edge, fraction 
of chord 

dimensionless sweep parameter, 


Vous modulus of elasticity, Ib/sq in. 

location of elastic axis measured from 
leading edge, fraction of chord 

dimensionless moment arm of section lift 
about elastic axis, e—a 

dimensionless moment arm of lift due to 
aileron deflection about elastic axis, 
срг— ё 

root-stiffness function given in equation 
(B25) of reference 2 

allowable bending stress, lb/sq in. 

dimensionless function of distance along 
span used in approximate formulas for 
angle of attack .due to aeroelastic 

' action of ailerons 

modulus of rigidity, Ib/sq in. 

wing thickness, in. 

section bending moment of inertia, in.* 

mass moment of inertia of entire airplane 
about its longitudinal axis, in.* 

mass moment of inertia of both wings 
about longitudinal axis of airplane, in.4 

section twisting moment of inertia, 1р.“ 

dimensionless parameters used in approx- 
imate formulas for dimensionless dy- 
namic pressures at aileron reversal 
given in table I 

dimensionless sweep parameter, 

2: E tan A; the -dimensionless 
ес, (ЕТ), 
parameter Ё/є is identical to k except 
that eis replaced by es 

lift per unit distance along span, lb/in. 

rolling moment on both wings, in-lb 

bending moment about ап axis perpen- 
dicular to elastic axis, in-lb 

free-stream Mach number 


N 


р 
pb[2V 
q 

q* 


Ri 


oo 


8* 


design load factor 

rolling acceleration, radians/sec? 
wing-tip helix angle due to roll, radians 
dynamic pressure, lb/sq ft 
dimensionless dynamic pressure, 


Cr, е1с;28 7 cos A | | | 
144 ШШ Т>; ES: Ше dimensionless 
А r 


dynamic pressure eg* is identical to 4% 
except that e; is replaced by e; 
dimensionless dynamic pressure, 
q Cua prsi sin А 
144 (ЕП, _ 
total wing area, sq in. 
distance along elastic axis measured from 
wing root, in. 
dimensionless distance along elastic axis, 
8/8: . 
accumulated torque about elastic axis, 
in-lb 
distributed torque due to inertia loading, 
in-lb/in. 
airspeed, ft/sec 
design gross weight of airplane, Ib 
weight of primary structure of both wings, 
Ib 
lateral coordinate, in. 
angle of attack in planes parallel to plane 
of symmetry, radians 
angle of attack equivalent to unit aileron 
deflection, radians | 
angle of local dihedral, radians; or span- 
wise slope of normal displacement of 
elastic axis 
aileron deflection measured in planes 
parallel to airstream, radians 
moment-arm ratio, 62/61 
structural-effectiveness factors defined in 
equation (15) of reference 2 
structural-effectiveness factors defined in 
table 1 of reference 2 
angle of sweepback at elastic axis 
wing taper ratio, с/с, 
density, slugs/cu ft 
angle of structural twist in planes pet- 
pendicular to elastic axis, radians 
tip stiffness ratio, an 
(E 'constant stress 


unit step function of distance along span 


at divergence 

at inboard end of aileron; inertia, in 
equation (36b) 

at outboard end of aileron; except in M, 

at aileron reversal 

at wing root 

structural (due to structural deformation) 

at wing tip 

rigid wing (for g*=e«q*=¢=0) 
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USE OF THE CHARTS AND APPROXIMATE FORMULAS 


SUMMARY OF METHOD AND SCOPE OF THE CALCULATIONS ON WHICH 
CHARTS AND APPROXIMATE FORMULAS ARE BASED 


А. brief description of the method and scope of the cal- 
culations is given here to indicate the limitations of the 
charts and approximate formulas. A detailed description 
of the method is given in the appendix. 

Most of these calculations were performed by an exten- 
sion, based on reference 1, of the matrix method described 
in appendix А of reference 2. "This method consists in 
solving the differential equations descriptive of an elas- 
tically deformed wing under aerodynamic loadings by 
numerical methods employing matrix techniques. "Treated 
by this method were wings with three taper ratios (1, 0.5, 
0.2), one aileron configuration (50 percent semispan, out- 
board), two types of stiffness distributions, several values 
of the sweep parameters Ё and d which include sweptfor- 
ward, unswept, and sweptback wings, and several values of 
the section moment-arm ratio e and the dynamic-pressure 
ratio q/qp. Calculated for each case were the dynamic 
pressure at aileron reversal and the changes in the spanwise 
lift distribution and rolling moment due to aileron deflection. 
For the constant-chord—constant-stiffness wings, calcula- 
tions were also performed by an extension of the analytic 
method described in appendix А of reference 2, which con- 
sists in solving the differential equations exactly for these 
relatively simple cases; these calculations were made for 
two aileron spans and several values of the parameters 
k, d, e, and 4/4р. In all cases the ratio of the aileron chord 
to the wing chord was constant along the span of the aileron. 

Some approximations have been made in the calculations 
concerning the aerodynamic induction effects, the root 
rotations, and the stiffness distribution, primarily in order 
to hold the number of variables considered in the analysis 
to а minimum and to make the results more generally 
applicable. 

Aerodynamic induction effects at subsonic speeds are 
taken into account by an overall reduction of the strip- 


theory loading and, in the matrix calculations, by rounding’ 


off the strip-theory loading at the tip (see refs. 1 and 2); 
for supersonic speeds strip theory is used with a small reduc- 
tion at the tip in the matrix calculations. This approxima- 
tion has made it unnecessary to consider explicitly the 
effects of aspect ratio, sweep, and Mach number on the 
rigid-wing lift distribution; the effects of these parameters 
on the total lift and on the aeroelastic increment to the lift 
distribution have been taken into account. 

The rigid-body rotations imparted to a swept wing by its 
triangular root portion vary among different designs in a 
largely unpredictable manner. They have therefore been 
taken into account only by the use of an effective root, the 
selection of which in any given case is discussed briefly in a 
subsequent section. 

The spanwise distributions of the bending and torsional 
stiffnesses depend on the detailed design of the wing and 
cannot be generalized easily. The stiffmess distributions 
used in the calculations of aeroelastic effects were obtained 
from the constant-stress concept outlined in appendix B of 
reference 2, which constitutes an effort to relate the stiff- 
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ness of a wing to ite strength on the basis of the following 
assumptions: 


(1) The combined bending and torsional stresses are 
constant along the span. 

(2) The bending and torsional stresses are combined in 
such a manner that the sum of the ratio of the actual to the 
allowable bending stress and the ratio of the actual to the 
allowable torsion stress is equal to unity when the margin 
of safety 1s zero. 

(3) The structure is of the thin-skin, stringer-reinforced 
shell type and its main features do not vary along the span; 
for instance, the number of spars and their chordwise loca- 
tions are constant along the span. 

(4) At the design condition the spanwise distribution of 
the applied loading is proportional to the chord. 

Also used in the calculations were stiffness distributions 
which vary as the fourth power of the chord, as do those of 
solid wings and wings with seometeically similar cross 
sections; as pointed out in a subsequent section, the results 
of these calculations can be used to estimate aeroelastic 
phenomena of some wings which have large cutouts or 
which for some other reason do not have stiffness distribu- 
tions represented fairly closely by those of the constant- 
stress type. 

All calculations are based on the ee that twist- 
ing is resisted primarily by the torsion cells of the wing 
structure and that the wing deformations can be estimated 
by means of the elementary theories of bending and torsion 
about an elastic axis. 

SELECTION OF PARAMETERS 
Geometric parameters.—The geometric parameters used 
in the analysis are defined in figure 1. The location of the 
Center of pressure due to aileron deflection... 
Elastic oxi. N 
Aerodynamic center. ` 


. Effective root 


-Elastic axis 


| NS Effective tip 
ғ с 
с 


Еісовв 1.—Definitions of geometric parameters. 


м 
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effective root indicated in this figure is discussed in refer- 
ence 2. In the present report the angle of aileron deflection 
6 is defined as being measured in planes parallel to the air- 
stream. This angle is equal to the product of the angle of 
rotation of the aileron about the hinge axis and the cosine of 
the sweep angle of the hinge axis. 

Although most of the charts апа approximate formulas 
are based on а half-span outboard aileron (8,*—0.5, 5,*— 1), 
the results of the analysis of the uniform wings with full- 
span ailerons (s,*—0, s,*—1) and “ыр ailerons” (8,*—1, 
8,*—1) indicate that, except for the angle-of-attack distri- 
butions, the results based on & half-span aileron may be 
expected to be valid for outboard ailerons having spans 


which differ considerably from one-half. 
Aerodynamic parameters.—The aerodynamic parameters 


which enter the analysis are the effective wing lift-curve 
slope, the location of the wing aerodynamic center, the loca- 
tion of the chordwise center of pressure due to aileron 
deflection, and the angle of attack equivalent to unit aileron 
deflection. An effective wing lift-curve slope Cz, , ар- 
plicable to basic lift distributions due to built-in twist, 
aileron deflection, roll, or aeroelastic twist, is approximately 
given at subsonic speeds by the relation 


A cos А 


C, e, —À (1) 
4c r45* cos A 
T 
where c;, is approximately given by 
2T (2) 


"a —Mj сов? A 
The basis of equation (1) is explained in reference 2. At 
supersonie speeds (more specifically, for supersonic leading 
and trailing edges), the effective wing lift-curve slope is 
approximately 
О, = 4 cos А 
ne JM? cos? A— 1 


provided М, is greater than 1/cos A. If M, is greater than 
1 but less than and not too close to 1/cos A, equations (1) and 
(2) may be used in the absence of better information ; however, 
the results obtained for this range of Mach numbers should 
be used with caution. 

The lift-curve slopes given by equations (1) and (3) should 
not be confused with the rigid-wing lift-curve slope or the 
damping-in-roll derivative; they are merely effective values 
suitable for aeroelastic calculations. Values of the rigid- 
wing lift-curve slope and the damping-in-roll derivative can 
be used in conjunction with the methods of the present 
report and of reference 2 to obtain the values of the flexible- 
wing lift-curve slope and the damping-in-roll derivative 
because means are presented herein and in reference 2 for 
estimating the ratios of the flexible-wing to the rigid-wing 
values. For this purpose any experimental information 
concerning the rigid-wing values can be used; if none is 


(3) 


available, references 3, 4, and 5 may be used at supersonic . 


speeds and reference 6, at subsonic speeds. 

The local aerodynamic centers are assumed to be at a 
constant fraction of the chord from the leading edge, so that 
they are all equal to the wing aerodynamic center as a fraction 
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of the mean aerodynamic chord. The moment arm e, із 
then given by the relation 


е=е-—а (4) 


The loeal'centers of pressure of the lift due to aileron de- 
flection are also assumed to be a constant fraction of the 
chord from the leading edge; and the moment arm e; is then 
given by 

€62— сре (b) 


Theoretical two-dimensional values of the parameter cp; are 
presented in figure 2 for both trailing-edge and leading-edge 
ailerons at subsonic and supersonic speeds. At subsonic 
speeds the effect of finite span is to shift the center of pressure 
rearward. An appropriate value for this rearward shift 
may be estimated from the following relation, based on 
lifting-line theory for unswept elliptic wings with full-span 
ailerons: 


Rear: 1 
Pig Pin (Png) 


where the subscripts IT and III refer, respectively, to two- 
and three-dimensional values. The use of the swept-span 
aspect ratio Ал in place of A should serve to extend this ap- 
proximate relation to swept wings. 

Theoretical two-dimensional values of оз, the angle of 
attack equivalent to unit aileron deflection, are also given 
in figure 2. At low aspect ratios the values of оз for sub- 
sonic speeds tend to be higher than these two-dimensional 
values; as the aspect ratio approaches 0, o; approaches 1, at 
least in the case of wings without reentrant trailing edges. 
Experimental values of both o; and cp, are preferable to 
theoretical values if they are available. For spoilers, ex- 
perimentally obtained values have to be used. 

The effective lift-curve slope and the values of the section 
moment arms vary with the free-stream Mach number; 
hence, the appropriate values must be used at each flight 
condition for which aeroelastic calculations are made. 

The airspeeds at which the lateral-control aeroelastic 


` phenomena are of interest enter the calculations in the form 


of the corresponding dynamic pressures. These dynamic 
pressures, in turn, are expressed in dimensionless form by 
means of the relations 

Сы, 6/82 сов А 


S Na EO ке 
177144 (GJ), (6) 

. C, £207 8? соз А 
«та (82) 9 

or 
Сі, ¢6 sin A 
gu Уа Ааваа (8) 
144 (ЕТ), 


The ratios of these quantities, 


EIN (GJ). tan A (9) 


je 
q ёс, (ЕТ), I) 


and 


LM кн EIN (GJ). 
Eu л Up s 


are independent of the dynamic pressure and are very useful 


(10) 
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for analyzing the aeroelastic behavior of swept wings. Two 
other dimensionless parameters, which are independent of 
the dynamic pressure, enter the calculations: 

(11) 


Е = — 


ё1 


(GJ). |. 4 
(ЕТ), m А 


Structural parameters.—For the purposes of an aeroelastic 
analysis the wing structure is characterized by the location 
of the elastic axis, the magnitude and distribution of the 
bending and torsional stiffnesses EJ and GJ, and the magni- 
tude:of the rigid-body rotations imparted to the wing by its 
root (taken into account in this report only by the location 
of an effective root). The selection of these structural 
parameters is discussed in reference 2. 

PRELIMINARY SURVEY OF LOSS IN LATERAL CONTROL 

The information contained in some of the charts and ap- 
proximate formulas presented in the following sections of 
this reporb has been summarized in figures 3 to 5 for the 
purpose of ascertaining in advance of more detailed esti- 
mates, if desired, whether the aeroelastic phenomena con- 


and 
(12) 


sidered in this report are likely to affect the design of the. 


wing structure. This preliminary survey is not essential 
to any of the further calculations but may show them to be 
unnecessary in some cases. "These figures pertain to constant- 
stress wings with half-span outboard ailerons. 

The charts of figures 3 (a), 4 (a), and 5 (a) pertain to wings 
of taper ratios 0.2, 0.5, and 1.0, respectively, with the moment 
arm ez equal to 0 (corresponding roughly to subsonic flow 


Ке 
о РЕ 
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ч е7 727 
Т 
НИЕ 

нн 

Е JE 

pas 

5 B Ши азы. ІШЕ 

8" LLL TL LEES 


4 


: 6 
Aileron -chord ratio, c; /c 
(а) Trailing-edge ailerons. 
Figure 2.—Theoretical two-dimensional values of the aileron-force parameters. 


821005--55----16 


reference 2. 


Aileron- effectiveness 


Center of pressure due to 
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conditions and elastic-axis locations fairly far back on the | 
wing). These figures show the dynamic-pressure parameter 
4% defined by either equation (6) or 
« 4 Cr, €i cos A { 
_ 18432 9 nW h, 
Fa S C т1а | 
plotted against the sweep parameter k defined by equation 
(9) or 


(13) 


Ал 
2 Е ш 71961 ип 
for several values of the aileron effectiveness parameter 
Cil Cry, and for two values of the stiffness parameter 
d _ (ЕГ), €1C; : 


1+6 


Ё (14) 


E (GJ. \ в, (5) 
or 
8а 892. Е лме!) 1617 (16) 
e yP G Ал 1879771515. 


where F, is a root-stiffness parameter, and та, 7», Ne, 77, 78, 
7$, 715, Dd 719 are structural-effectiveness factors defined in 
(If, at the time a preliminary survey of аего- 
elastic effects is to be made, no information whatever con- 
cerning the wing stiffness is available, eqs. (13), (14), and 
(16) may be used; if an estimate of the root stiffnesses 
(GJ), and (EI), is available, eqs. (6), (9), and (15) may be 
used.) 

The charts of figures 3 (b), 4 (b), and 5(b) are the same as 
those of figures 3 (а), 4 (а), апа 5 (а), respectively, except 
that these charts are for wings with the moment-arm ratio 
e equal to unity (corresponding roughly to subsonic conditions 
with the elastic axis fairly far forward on the wing). 
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(b) Leading-edge ailerons. 


Oynomic-pressure parameter, 9* 


Dynamic-pressure роготвіег, g* 
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FrcunE 3=—Charts for a preliminary survey of lateral control for flexible wings of taper ratio 0. 2. 
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(b) Wings with moment-arm ratio в = 1.0. 
Бістнв 3.—Continued. 
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(c) Wings with moment arm e,—0. 
ЕтзовЕе 3.—Concluded. 
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(а) Wings with moment arm е=0. 
FraunE 4.—Charta for a preliminary survey of lateral control for flexible wings of taper ratio 0.5. 
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The charts of figures 3 (c), 4 (c), and 5 (c) pertain to wings 
with the moment arm є; equal to zero (corresponding roughly 
to supersonic flow conditions). These figures show the 
dynamic-pressure parameter eg* defined by either equation 
(7) or 

(4X? Cr, Lp cos A 
= 18432 @ пў м 
Fs 75. га 


plotted against the sweep parameter k/e defined by equation 
(10) or 


k 1-Х а Ал 
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for two values of the stiffness parameter 
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The various lines of the charts of figures 3 to 5 designate 
the conditions at which a wing designed on the basis of 
strength considerations alone is likely to encounter changes in 
aileron rolling moment by various amounts due to wing 
flexibility. The lines for zero rolling moment designate 
the conditions at aileron reversal. These charts should be 
used in conjunction with the preliminary survey charts in 
reference 2. The significance of the four quadrants of 
figures 3 (a), 3 (b), 4 (a), 4 (b), 5 (a), and 5 (b) is the same 
as that of the four quadrants of figures 2 (a), 2 (b), 2 (c), 
and 2 (d) of reference 2 and is discussed in the section “Рге- 
liminary Survey of Aeroelastic Behavior" of reference 2. 
The significance of the quadrants of figures 3 (с), 4 (с), 
and 5 (c) of the present report can be analyzed in the same 
way, except that the moment arm e; takes the place of the 
moment е. These three figures are analogous to figure 2 (e) 
of reference 2 in that they pertain to the case e;—0. (Al- 
though Ё is infinite when e, is zero, which happens when the 
section aerodynamic centers are on the elastic axis, the 
parameter k/e is not infinite in that case, except if e; also 
happens to be zero, a condition which can be realized only 
with a full-chord aileron.) 
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After № has been found, through the use of the charts in 
reference 2, to what extent the wing design is affected by 
such aeroelastic phenomena as divergence and aerodynamic- 
center shift, the same procedure can be used with the pre- 
liminary survey charts in this report to ascertain whether 
the wing design is likely to be affected by lateral-control 
difficulties. If these charts indicate the likelihood of sig- 
nifieant aeroelastic effects on the aileron rolling moment, 
further calculations are desirable. The charts and approxi- 
mate formulas of this report may be used for the preliminary 
calculations; once the structure has been designed, more 
refined methods, such as that of reference 1, may be used. 
CALCULATION OF THE AEROELASTIC PHENOMENA RELATED TO LATERAL 

CONTROL | 

Analysis of the many solutions for the aeroelastic phe- 
nomena considered in this report obtained by the methods 
given in the appendix shows that the data can be sum- 
marized by means of approximate formulas. These for- 
mulas involve the aerodynamic, geometric, and structural 
parameters of the wing through the dimensionless param- 
eters С, k/e, e, and d (eqs. (9), (10), (11), and (12), respec- 
tively) and through a series of constants К, to К; The 
constants are functions of the taper ratio, stiffness distribu- 
tion, and aileron span and are given in table I. As in 
reference 2, the form of these approximate formulas has 
been guided by considerations based on an idealized semi- 
rigid wing, and the actual values of the constants K, to 
К, (К, and К» have been given in ref. 2) were obtained by 
fitting the solution for the functions В, to В, defined in the 
appendix by equations (A18) and (A27) to their approximate 
expressions, equations (A32). 

Dynamic pressure at aileron reversal.—The solutions for 
the aileron reversal speed obtained by the methods given in 
the appendix can be summarized by approximate formulas 
which give the dimensionless parameters q*r, (eg*)s, or 
4к--іһаб is, the values of the parameters defined in equa- 
tions (6), (7), and (8) which correspond to the value of the 
dynamic pressure at aileron reversal—in terms of the 
parameters Ё, k/e, e, and d defined by equations (9), (10), 
(11), and (12), respectively. 

An approximate formula for @*в is 
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For very small values of the section moment arm e, and 
the resulting large values of the parameters e and k, the 
following alternative forms of equation (21) are more 
convenient to use: 


к (1—к.; а) 
Т. (22) 
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When the angle of sweep is zero, equations (21) and (22) 
reduce, respectively, to | 


(24) 


(25) 


and when the moment arm е, is zero, as it may be in super- 
sonic flow, equations (22) and (23) reduce to’ 
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The constants K, to K, are given in table I for wings 
having half-span outboard ailerons and taper ratios of 0.2, 
0.5, and 1.0, for the two different types of stiffness distribu- 
tions. These constants were found from the results of the 
numerical matrix method derived in the appendix. Also 
given in table I are these constants for uniform wings having 
full-span ailerons and tip ailerons calculated from the results 
of the analytic integration method of the appendix. Since 
values of the constants K, to К, are given for three aileron 
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spans in the case of the uniform wing, they may be inter- 
polated to yield values for other aileron spans. Мо calcu- 
lations were made for other than half-span ailerons on 
tapered wings; nevertheless, as pointed out previously, qa 
calculated for half-span ailerons should be reasonably valid 
for outboard ailerons having spans which differ considerably 
from one-half. 

No aileron-reversal calculations have been made for swept 
wings with inboard ailerons. However, for unswept uni- 
form wings the dynamic pressure at aileron reversal has been 
calculated for the limiting case of a wing with ап mboard 
aileron of vanishingly small span by operational methods 
similar to those described in the appendix in connection 
with the calculations for uniform wings. The value of g*z 
obtained in this manner is shown as a function of e in figure 
6. Also shown in figure 6 are the values of 4% for full-span 
and for tip ailerons. For small values of e, such as are likely 
to be encountered in subsonic flight, there is little difference 
in the values of 4% for the three aileron configurations, but 
at large values of e, such as are likely to be encountered at 
high-supersonic speeds, there is some difference between 
them; the aileron reversal speed is highest for the tip aileron 
and lowest for the inboard aileron. However, these con- 
clusions may not be valid for nonuniform or swept wings. 

With the values of g*s, (eg*)r, and qz given by equations 
(21), (22), апа (23) and the definitions of these parameters 
given by equations (6), (7), and (8), the values of q at 
aileron reversal may be determined. If desired, the corre- 
sponding airspeed may be determined from the relation 


=-/-® 
Va р[2 on 


If both gp (as obtained from ref. 2, for instance) and gz 
are positive and дв is greater than qp, there 18 no actual 
aileron reversal speed; in fact, the rolling moment due to 
aileron deflection will tend to increase with dynamic pressure 

until divergence is reached. 
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Figure 6.—Comparison of dynamic pressures at aileron reversal 
calculated by the analytic integration method of the appendix 
with those calculated from equation (A36) for unswept uniform 
Wings. 
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Nor is there an actual reversal speed if gp is negative, 
regardless of whether (р is positive or negative, although the 
rolling moment due to aileron deflection will decrease slightly 
with increasing dynamic pressure if gg/qp 18 greater than 1. 
Inasmuch as aileron reversal, unlike divergence, is not an 
instability problem, although it is often convenient to analyze 
it as such, there is no possibility of encountering aileron 
reversal in a higher mode when qz is negative, at least, not 
in the case of ordinary wings with straight leading and 
trailing edges and with substantially straight elastic axes. 

The value of gp calculated for any given value of q*z, 
(eq *)5, ог Gr depends on the value of the effective lift-curve 
slope Cra, and, hence, on the Mach number. As suggested 


in reference 1, the value of qr may be plotted against Mach 
number on log-log coordinates; if the straight lines of the 
actual dynamic pressure at several altitudes as functions of 
Mach number are drawn on the same plot, an intersection 
of the reversal line with one of the lines of actual dynamic 
pressure designates possible aileron reversal at that value of 
dynamic pressure, Mach number, and altitude. | 

Spanwise angle-of-attack distribution.—In the appendix, 
an approximate expression is determined for the change in 
angle of attack due to the deflection of ailerons on flexible 
wings. The ratio of the angle-of-attack distribituon due 
to structural deformation o, to the effective angle of attack 
of an aileron озб 18 


саб 1 _ 9, 1— К, 
dp 


Тһе functions f, and Af, of the spanwise coordinate 8% are 


given in figure 7 for wings having half-semispan outboard 
ailerons, taper ratios of 0.2, 0.5, and 1.0, and the two different 
typés of stiffness distribution. The value of gp required in 
equation (29) may be found from the approximate expression 
for q*» or 9р given in reference 2 as 


‚+ _ К 
or 
o=- (81) 
Е 


When др is very large, а more convenient form of equation 
(29) is as follows: 


“ _ (32) 
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— 


dp 
Spanwise lift distribution Within the limitations of the 
modified strip theory used in the analysis, the lift per inch 
of span is proportional to the local effective angle of attack, 
so that 


l а; 
а. 4601, (6+1) (38) 
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where a,/as5 is obtained as indicated in the preceding section 
and I, is а unit step function of the distance along the span 
defined by 
1,70 (when 8< 81) 
1,=1 (when s>s,) 

Inasmuch as the functions f, and Af, have been calculated 
only for half-semispan outboard ailerons or spoilers, the 
expressions for angle-of-attack distributions and lift. distri- 
butions (eqs. (29) and (83)) can be used directly only for 
this case. 

Rolling-moment coefficient and rate of roll due to aileron 
deflection.—The rolling-moment coefficient due to a unit 
aileron deflection C;, may be obtained in terms of its equiva- 
lent rigid-wing value from the approximate formula 


q 
Ci dr 
Zi (84) 
Cry, 1-1. 
9р 


апа Ше wing-tip helix angle due to an aileron defection 
pb/2V, which is a measure of the rate of roll and the rolling 
maneuverability, may be obtained from 
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(b) Stiffnesses related to those given by constant-stress criterion for 
wings of taper ratio 0.5. 
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The manner in which the aileron-effectiveness parameter 
Cl C varies with dynamic pressure depends on the ratio 
of the reversal to the divergence dynamic pressure as may 
be seen from equation (34) and figure 8. When the aero- 
dynamic center is ahead of the elastic axis, as is generally 
the case at subsonic speeds, q&/gp is positive and greater 
than one for sweptforward wings, positive and less than one 
for unswept wings, and negative for sweptback wings. 
Thus, in general, бб, increases with dynamic pressure 


until divergence is reached for sweptforward wings; it de- 
creases slowly at first and then more rapidly as reversal is 
approached for unswept wings; and it decreases rapidly at 
first and then more slowly for sweptback wings. The rapid 
decrease in rolling effectiveness for sweptback wings can be 
alleviated by the use of unconventional lateral-control 
devices which have their centers of pressure ahead of the 
elastic axis, such as leading-edge ailerons or spoilers. These 
devices may also serve to make the dynamic pressure at 
aileron reversal negative, so that there is no reversal of 
lateral control in the given speed range. The loss or gain 
of lateral control is then given by the part of figure 8 for 
negative values of 4/4в. 

The analysis summarized by the approximate formulas 
(34) and (35) is based on rolling moments about the wing 
root instead of about the fuselage center line, partly to 
simplify the analysis and partly to avoid the introduction of 
the fuselage width as another independent parameter. How- 
ever, these approximate formulas should be valid for ob- 
taining rolling moments and rates of roll about the fuselage 
center line as well, because equations (34) and (35) are ex- 
pressed as ratios of flexible-wing to rigid-wing values; that 
is, the ratio of the rolling moment about the fuselage center 
line to its rigid-wing value should be nearly the same as the 
ratio of the equivalent rolling moment about the wing root 
to its rigid-wing value when the ratio of the fuselage width 
to the wing span is small. 

The rolling-moment coefficient and rates of roll given by 
equations (34) and (35) are functions of the aileron span 
inasmuch аз gz is a function of the aileron span. The varia- 
tion of C;, or pb/2V with aileron span can therefore be found 
only for uniform wings; however, since the effect of aileron 
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Fiaurp 8.— Variation of aileron rolling effectiveness with dynamic- 
pressure ratio and the parameter qr/gp. 
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span on gr is Dot very great, its effect on С, or pb/2V is not 
likely to be great. The rolling moment due to the defloction 
of an mboard aileron can be found by superposition, because 
the aeroelastic equations upon which the results of this ro- 
port are based are linear (that is, the rolling-moment coeffi- 
cient due to a 30-percent-span inboard aileron, for example, 
is equal to (1, for a full-span aileron minus Сі, for a 70- 
percent-span outboard aileron), but only in the case of the 
uniform wing is the required information for the full-span 
aileron presented herein. That the aeroelastic effects on 
lateral control are likely to be similar for inboard as for out- 
board ailerons, at least for unswept wings, may be deduced 
from figure 6. 

Inertia effects.—In steady rolling flight no inertia- effects 
are present which can affect the static aeroelastic problem 
except, possibly, for centrifugal forces on heavy underslung 
nacelles. The maximum value of pb/2V is therefore usually 
unáffected by inertia effects. However, in the equally im- 
portant problem of initial rolling acceleration, which governs 
the time in which & given rolling velocity can be attained, 
inertia effects must usually be taken into account. In 
reference 2 the observation was made that in symmetric 
flight inertia effects are not as important as other static 
aeroelastic effects, except for flying wings, because the inertia 
forces are in about the same ratio to the aerodynamic forces 
as the wing weight is to the airplane weight. By the same 
reasoning, however, inertia effects are almost always very 


.important in getting into a roll, because the inertia forces 


are then in about the same ratio to the aerodynamic forces 
as the moment of inertia of the wings about the longitudinal 
axis of the airplane is to the moment of inertia of the entire 
airplane about its longitudinal axis, a ratio which is usually 
not much less.than 1. ' 

No charts are presented in this report for these inertia 
effects because the manner in which mass is distributed 
varies so widely among different airplanes that preparation 
of a generally applicable set of charts appears to be im- 
practical at present. However, the procedure outlined in 
reference 2 for taking mertia effects into account in the cal- 
culation of quasi-static aeroelastic phenomena by means of 
the charts presented therein may be applied to the calcula- 
tion of the inertia effects encountered in starting a roll. 
This procedure is described in the following paragraphs. 

For a given rolling acceleration р, the linear normal accel- 
eration of an element of mass at а distance y from the center 
line of the airplane is py. From this linear acceleration and 
the known or estimated mass distribution of the wing, the 


inertia load l; per inch of span and the inertia torque % per 


inch of span can be calculated for any given normal, pitch- 
ing, or rolling acceleration. Substitution of these loads and 
torques for the terms / and lec in equations (АЗ) or (А86) 
and equations (А2) or (A35) of reference 2, respectively, 
yields the values of the accumulated bending moments and 
torques in equations (A4) and (A5) or in equations (À37) 
and (A38) of reference 2. Equation (Аб) of reference 2, or 
the matrix equivalent of this equation, then yields the angle- 
of-attack distribution due to the deformations caused by the 
inertia effects associated with the given acceleration. 
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This angle-of-attack distribution can be considered as а 
geometrical angle-of-attack distribution; for the purpose of 
calculating the increment caused by aeroelastic action, this 
distribution can be approximated by a linear-twist angle-of- 
attack distribution with а value &t the wing tip which is such 
that the moment about the effective wing root of the area 
under the linear-twist distribution equals the moment of the 
area under the calculated angle-of-attack distribution due to 
inertia effects. (The moment, rather than the area, is 
suggested as а basis of correlation because the angles of 
attack near the wing tip are more important in aeroelastic 
phenomena than those at the wing root.) The justification 
for this rather arbitrary approximation to the angle-of- 
‘attack distribution is that the correction to be applied as a 
result of aeroelastic action to the deformations due to inertia 
loads is usually small compared with these deformations. 

The angle of attack due to structural deformation a, 
associated with the linear-twist distribution can then be 
obtained from equation (21) and figure 7 of reference 2 or, 
if \=0, from figure 8 of reference 2. The lift distribution 
associated with the total angle-of-attack distribution due 
to the deformations caused by the inertia effects, including 
the increment in this angle-of-attack distribution produced 
by aeroelastic action, can then be found from equation 
(24b) of reference 2, in which о, and lh pertain to the calcu- 
lated angle-of-attack distribution due to the inertia effects 
(not the linear approximation to this distribution). This 
lift distribution ean be integrated to obtain the rolling 
moment due to inertia effects, as modified by aeroelastic 
action. | 

The rolling moment calculated in this manner may then 
be combined with the rolling moment due to aileron or 
spoiler deflection, which шау be calculated as indicated in 
. the preceding section. If the contributions of the tail and 
the fuselage to the rolling moment are neglected, 
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where f is the angular acceleration in roll, 7 is the mass 
moment of inertia of the entire airplane about its longitudinal 
axis, Т„ is the mass moment of inertia of both wings about 
the longitudinal axis of the airplane, and С, is the flexible- 


wing value of the rolling-moment coefficient due to aileron 
deflection (which may be calculated in the manner de- 
scribed in the preceding section). The ratio (OL',/Op), is 
the roling moment per unit rolling acceleration due to 
inertia effects, including aeroelastic effects, and is equal to 
—1, plus the rolling moment due to the lift which results 
from the deformations due to, the inertia loads per unit 
angular acceleration in roll as well as from the aeroelastic 
deformations which accompany these inertia deforma- 
tions; in other words, (OL’,,/07), is equal to —J,, plus the 
roling moment calculated as described in the preceding 
paragraphs for p—1. Then 


C;, 69 S b 
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where | 
1 
Crs, , E 1 OL’ Crs, 
І---- =) 
I—I,,\ ФР /s 
is a rolling-moment coefficient per unit aileron deflection 


which includes static aeroelastic effects, inertia effects, and 
the a&eroelastic magaification of the inertia effects. This 
rolhng-moment coefficient is a truer index of the rate at 
which a roll can be initiated than Са The initial rolling 


acceleration (disregarding  uüsteady-lift effects) can be 
calculated from equations (36). 


ILLUSTRATIVE EXAMPLE 


The approximate formulas described in the preceding 
sections have been used to find the effects of aeroelasticity 
on some lateral-control properties of the wing considered in 
the illustrative example of reference 2. The resulting cal- 
culations are ап extension of those in reference 2, and the 
&dditional parameters are presented in table II. 


TABLE IIL.—PARAMETERS OF EXAMPLE WING 


Bu ріс 
(М=з1.5) 


Subsanio (2-0) 
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The subsonic and supersonic values of the parameters Ё, 
є, and @ were calculated from equations (9), (11), and (12), 
respectively. With the appropriate values of the factors K, 
to К, interpolated from table I, the values of @*ь were cal- 
culated from equation (21) and are given in table IT. From 
these values of 4%;, the subsonic and supersonic dynamic 
pressures at aileron reversal were found by means of equation 
(6). These values of gg vary as the reciprocal of the effective 
lift-curve slope, if the corresponding values of e; and e; are 
assumed to remain constant. 

In order to find the angle-of-attack distributions due to 
deflections of the aileron from equation (29), the values of 
the functions ў, and Af, were taken from figure 7(c). The 
spanwise change in angle of attack is shown in the top plot 
of figure 9 for different values of the dynamic-pressure ratio. 
The rolling-moment coefficient and wing-tip helix angle due 
to deflections of the aileron were calculated from equations 
(34) and (35) and were plotted in figure 9 as functions of the 


dynamic-pressure ratio 9 . і 
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Етаунв 9.—Effect of aeroelastic action on some lateral-control prop- 
erties of the example wing. 
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DISCUSSION 
LIMITATIONS OF THE CHARTS AND APPROXIMATE FORMULAS 


The charts and the approximate formulas presented in 
this report are subject to certain limitations as a result of 
the approximations made in the calculations on which they 
are based. "These limitations are discussed fully m reference 
2 and can be classified as restrictions on the plan form, on 
the speed regime, and on the wing structure. The limitations 
are given very briefly as follows: 

(1) The results obtainable by the use of the charts and 
approximate formulas are likely to be unsatisfactory for 
wings of very low aspect ratio, very large sweep, or zero 
taper ratio. 

(2) The results are restricted to wings on which the span- 
wise lift distribution is roughly proportional to the chord 
and angle of attack and on which the section acrodynamic 
centers are at an approximately constant fraction of the 
chord; these restrictions are most likely to be violated by 
wings flying at transonic speeds and by wings having con- 
centrated sources of lift, such as nacelles and tip tanks. 

(3) The results are somewhat restricted to wings with one 
of the two types of spanwise stiffness distributions used 
in the analysis, wings with no chordwise bending (relatively 
thick wings), and wings having an elastic axis at an approxi- 
mately constant fraction of the chord. 

The manner in which the aeroelastic effects of aileron 
deflection are analyzed in this report imposes certain ad- 
ditional limitations that particularly affect the aileron 
geometry. In the analysis in the appendix, the spanwise 
lift distribution due to the deflection of an aileron was 
approximated by strip theory, an approximation which is 
probably less valid for aileron deflections than for geometric 
angles of attack. The assumption was also made that the 
centers of pressure due to aileron deflections are at a constant 
fraction of the wing chord; this assumption is also probably 
less valid than the assumption that the section aerodynamic 
centers of the wing are at a constant fraction of the chord. 


Since these assumptions are more nearly true for wings of 


high than those of low aspect ratio, these limitations serve, 
in effect, to restrict the applicability of the present report 
to aspect ratios somewhat higher than those amenable to 
the analyses of reference 2. 

The results of the present report do not take into account 
explicitly any flexibility of the aileron itself because of the 
assumption that the angle between the aileron and wing is 
constant along the span of the aileron. This assumption is 
almost universally made in analyzing the aeroelastic prop- 
erties of ailerons and is justifiable because the net effect of 
the difference between the wing deformations and the aileron 
deformations on the overall lift and moments appears to 
be négligible. 

As a result of the fact that the static aeroelastic phenomena 
associated with lateral control involve many more param- 
eters than do those associated with symmetric flight, the 
coverage of the various parameters is not as complete as 
in reference 2. Specifically, all the charts and approximate 
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formulas presented herein are restricted to outboard lateral- 
control devices (ailerons or spoilers) except for the uniform- 
wing case, and most of the calculations upon which these 
results are based were made for wings with half-semispan 
outboard ailerons (8;*=0.5, s,*—1). However, the results 
of the analysis for the uniform wings with full-span ailerons 
($;,*=0, 8,*=1) and tip ailerons (s,*—1, s,*=1) indicate 
that, except for the angle-of-attack distributions, the results 
based on a half-span aileron are approximately valid for 
outboard ailerons of spans differing considerably from one- 
half. The ratio of the aileron chord to the wing chord is 
assumed in the charts and formulas to be constant over the 
span of the aileron. 


RELATION BETWEEN STRENGTH AND STIFFNESS AS DESIGN CRITERIA 


The relation between strength and stiffness as design 
criteria was discussed in reference 2. The preliminary sur- 
vey charts in reference 2 indicate the extent to which a 
wing designed on the basis of strength considerations alone 
is likely to be affected by aeroelastic phenomena and, con- 
sequently, indicate whether the wing has to be stiffened 
beyond the amount associated with the required strength. 
The preliminary survey charts of this report (figs. 3 to 5) 
serve the same purpose in regard to the aeroelastic effects on 
lateral control; furthermore, even though the charts of 
reference 2 may indicate that a particular wing is not sig- 
nificantly affected by the aeroelastic phenomena considered 
in that report, this wing may still have to be stiffened be- 
cause of an undesirably large loss in lateral control. 

As may be concluded from the survey charts of reference 2 
and the present report, as well as from the discussion: con- 
tained in reference 2, the following wings designed on the 
basis of strength considerations are most likely to be subject 


to adverse aeroelastic effects on lateral control and rolling. 


maneuverability: 

(1) Wings operating at a high speed or dynamic pressure 

(2) Sweptback wings 

(3) Thin wings 

(4) Wings designed for low wing s loading 

(5) Unswept and moderately swept wings with an elastic 
axis relatively far forward on the chord or with the center 
of pressure of the lift produced by aileron deflections rela- 
tively far back on the chord as a result of the aileron configu- 
ration (small aileron chord or wing of low aspect ratio) or 
flight condition (supersonic speeds) 

(6) Wings with a relatively high lift-curve slope 


STRUCTURAL WEIGHT ASSOCIATED WITH THE REQUIRED STIFFNESS 


When a given wing has been shown to be subject to un- 
desirably large aeroelastic effects by means of the charts of 
reference 2 and this report or.by any other method, the 
problem arises how to distribute the additional required 
stiffness—that is, which spanwise distribution of structural 
material will alleviate the adverse aeroelastic phenomena 

to the desired extent with the minimum increase in struc- 
tural weight. 

In order to shed some light on this problem, aeroelastic 
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and weight calculations have been made for wings with а 
taper ratio of 0.5 with a family of somewhat arbitrarily 
selected stiffness distributions which differ from the distri- 
bution required by the constant-stress criterion in a manner 
described in reference 2. These stiffness distributions are 
designated by the tip stiffness ratio w, which is the ratio of 
the stiffness EI or GJ at the wing tip to the corresponding 
stiffness of a constant-stress wing. ‘The results of the 
lateral-control calculations for wings with taper ratio 0.5, 
with constant wing-thickness ratio М/с along the span, and 
with two of these stiffness distributions are included in 
table I and figure 7 (b). The designation ‘‘excess strength" 
refers to the stiffness distribution increased over the constant- 
stress requirement to such an extent that the value of w is 
2.0. The results of the aeroelastic calculations for the 
stiffness distributions decreased below the constant-stress 
requirement to a value of w=0.5 happen to be the same as 
the results for the constant-stress stiffness distributions for 
wings with linearly varying wing-thickness ratio and 
ұлы --0.5. The structural weight considered in these cal- 
culations is that of the primary load-carrying structure; 
the remaining structure is assumed to be unchanged in the 


stiffening process. 
The results of the weight calculations and aeroelastic 


calculations in reference 2 indicated that in the case of wings 
with taper ratio 0.5 stiffening the structure in the outboard 
region of the wing (w greater than 1) was more efficient, 
from weight considerations, in alleviating the aeroelastic 
effects considered in that report than the addition of stiffness 
in the inboard region. This conclusion is corroborated, in 
essence, by the calculations made for the aeroelastic phe- 
nomena, considered herein. Figure 10, which consists of a 
plot of the structural weights required for a given loss (20 
percent) in lateral control at a given dynamic pressure, indi- 
cates that the least weight is associated with values of the 
tip stiffness ratio w greater than 1, except for wings with 
values of the sweep parameters Ё or k/e equal to —8. These 
large negative values of k or k/e, however, pertain. to wings 
that are either (1) sweptforward, or (2) sweptback with the 
aerodynamic center behind the elastic axis (in the case of 
negative Ё), or (8) sweptback with the center of pressure due 
to aileron deflection ahead of the elastic axis (in the case of 
negative k/e, as it may be for spoilers or leading-edge 
ailerons). 

For sweptforward wings, the aileron rolling moment 
usually increases rather than decreases with dynamic pres- 
sure, so that lateral control does not impose any structural 
requirements. For sweptback wings with negative values of 
e; and еҙ, alleviation of the aeroelastic effects in lateral con- 
trol can be effected at the least cost in weight by adding 
structural material in the inboard region of the wing; inas- 
much as the minimization of the shift of the aerodynamic 
center of. these wings can be effected most efficiently by 
stiffening the outer region of the wing, a compromise must 
be made if both types of static aeroelastic phenomena are of 
concern. 
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FiaurE 10.— The effect of tip stiffness ratio on the structural weight | 


required to maintain a constant level of lateral-control effectiveness 
(Cn —0.8C4,) at a given dynamic pressure for wings with a taper 


ratio of 0.5. 
SOME REMARKS CONCERNING THE AEROISOCLINIC WING 


As shown in reference 2 an overall type of aeroisoclinicism, 
in which bending and torsion action tend to cancel for the 
wing as a whole, can be achieved for the aeroelastic phe- 
nomena considered in that report by a choice of a suitable 
ratio of the bending to the torsion stiffness or by a choice of 
the elastic-axis location, that is, by satisfying the relation 


в. (GJ). 


ес; (ET), 
(where К, is given in table Y). However, also pointed out 


1 
{ап A= (37) 


REPORT 1139—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


in reference 2 was the fact that even if the conditions of 
equation (37) are achieved there may still be great losses in 
lateral control, and the wing may still be subject to adverso 
dynamic phenomena; in fact, the severity of adverse aero- 
elastic effects in the lateral control and of certain dynamic 
phenomena may be increased as a result of achieviog аего- 
isoclinicism. i | 
The results of the present report corroborate the conclu- 
sion concerning aeroelastic effects on lateral control. How- 
ever, by suitable additional modifications a wing which has 
been made aeroisoclinie can also be made to suffer no loss in 
lateral control due to aeroelastic action. As may be seen 
from figures 3 to 5 or from equation (34), the condition for 
no losses in lateral control is that дь be equal to qp or, by 
setting equation (21) equal to equation (30), that 


Ке (Ks-- Ez) E 4- Kd -(G-- Ку) т d=0 (88) 
' of which the condition 
Be AU Dr oeil да. (39) 


баст (ЕГ), К; 


usually is an approximate solution. "Therefore, in order to 
satisfy both of the conditions specified in equations (37) and 
(38), the section moment arm е, must be nearly equal to 
—ё or, in other words, the center-of-pressure parameters а 
апа cp, must be nearly equal, а condition which can be 
satisfied by using full-chord ailerons (all-movable wing tips) 
or а combination of geared leading- and trailing-edge ailerons, 
for instance. However, as pointed out in reference 2, 
attempts at solving static aeroelastic problems by aiming at 
aeroisoclinicism may tend to aggravate certain dynamic 
phenomena. The same statement must also be made con- 
cerning the foregoing methods of alleviating static aeroelastio 
effects on lateral control; these methods may, for instance, 
lead to flutter difficulties which may require excessive mass 
balancing. 


RELATION OF THE CHARTS TO DESIGN PROCEDURE 


The conventional procedure of designing a wing on the 
basis of strength requirements and later checking it for 
aeroelastic effects can be facilitated at several stages by 
using the methods described in the present report and in 
reference 2. Аз pointed out in reference 2, for instance, the 
preliminary-survey charts presented therein can be used to 
establish some static aeroelastic characteristics that would 
be obtained in symmetric flight if the wing were designed for 
strength alone. 

Tf these characteristics are deemed satisfactory, the de- 
sign can proceed on the basis of strength requirements alone. 
Tf, on the other hand, they are considered unsatisfactory, 
the wing must be stiffened. The amount of - additional 
material required can be estimated, as indicated in reference 
2, by interpolating between the results presented therein 
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for the constant-strength case, the “excess strength" саве, 
(Ме), — 
(he), 
additional structural material is usually most effective if 
distributed near the wing tip. 

Similarly, the preliminary-survey charts of this report 
can be used to ascertain whether the wing can meet lateral- 
control requirements if designed for strength alone. If it 
must be stiffened to meet these requirements, the necessary 
amount of additional material can be estimated in the 
same manner as indicated in reference 2 for aeroelastic 
effects incurred in symmetric flight. 

Inasmuch as the charts of reference 2 and of the present 
report pertain only to static aeroelastic phenomena, the 
problem remains of ascertainiog in the preliminary design 
stage whether a wing designed for strength alone (or, for 
that matter, а wing designed both on the basis of strength 
requirements and of static aeroelastic considerations) is 
likely to experience flutter difficulties. However, flutter is 
а much more complicated phenomenon and depends on 
many more parameters than do static aeroelastic phenomena. 
Consequently, preparation of & generally applicable set of 
charts appears’ impractical. Nonetheless, although phe- 
nomenologically flutter is not related to the static aero- 
elastic phenomena considered in reference 2 and herein, it is 
mechanically related by virtue of the fact that all these 
phenomena depend on the wing geometry and the wing 
stiffness (although the aerodynamic parameters are different 
and flutter, unlike the static aeroelastic phenomena, involves 
the mass distribution of the wing and the damping properties 
of the structure). On the basis of past experience, certain 
qualitative conclusions can be drawn concerning this 
relation. 

As shown in the charts of this report, the aileron reversal 
speed, or the speed at which a specified amount of control is 
retained, is lower for highly sweptback wings than for 
unswept wings. Similarly, the divergence speed decreases 
rapidly as the angle of sweepforward increases. For a 
typical family of wings the values of g* at divergence and 
at reversal as obtained by the charts herein are shown as a 
function of the sweep parameter Ё in figure 11. For un- 
swept wings the dynamic pressure at flutter is usually within 
a certain range varying between a value lower than the 
dynamic pressure at reversal to a value higher than the 
dynamic pressure at divergence, depending on the geo- 
metric, structural, aerodynamic, and mass parameters of 
з given case, and varying even for a given case and a 
given speed range with altitude, because a change in air 
density may change the mode in which the wing flutters. 
This wide range is indicated in figure 11 by starting three 
flutter curves (which do not necessarily describe the upper 
and lower limits) at k=0. If the variation of flutter speed 
with sweep angle is assumed for the purpose of illustration 
to be similar to that indicated in figure 17 of reference 7, 


and the case 


As previously mentioned, the 
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the flutter curves of figure 11 are obtained. 

This figure must not be construed as presenting any 
quantitative information; to emphasize this point the family 
of wings is not identified. Even qualitatively the relation 
between the dynamic pressures at flutter, divergence, and 
reversal is subject to certain limitations because the flutter 
tests of reference 7 were-performed.at subsonic -speeds on 
models without ailerons and concentrated masses, which 
fluttered in the classical two-degree-of-freedom mode. There 
is reason to believe that sweptback wings with high aspect 
ratios flying at high altitudes may ‘experience a possibly 
mild form of flutter in a single-degree-of-freedom mode, be- 
cause a.vertical motion necessarily implies vertical bending 
and, hence, in the case of a swept wing, а variation in the 
angle of attack. In general, the greater the number of 
degrees of freedom the more difficult it is to relate flutter to 
the static aeroelastic phenomena. 

However, at subsonic speeds/and low or moderately high 
altitudes at least, the trend shown in figure 11 should be valid 
for wings without very large concentrated masses and with 
irreversible controls, which tend tò minimize the possibility 
of aileron-coupled flutter. Consequently, if these wings are 
highly swept back they can be designed to meet lateral- 
‚ control requirements with the likelihood that they will then 
be safe against flutter as well, provided conventional lateral- 
control devices are used. On the other hand, if these wings 
are unswept or even moderately swept back, they may have 
to be stiffened beyond the amount required by static aero- 
elastic considerations, or mass balanced, or both. 

In any event the final design must be checked both for 
static aeroelastic effects and for flutter. In many cases the 
static aeroelastic effects can probably be calculated with 
sufficient accuracy by means of the charts and approximate 
formulas of reference 2 and the present report. In some 
cases, however, particularly if these effects are in any way 
eritical, & more refined method of analysis, such ss that of 
references 1 and 8, may have to be used. 
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Етатвв ll.—Variation of critical dynamic pressures with the sweep 
parameter k. 
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CONCLUDING REMARKS 


An approximate method based on charts and approximate 
formulas has been presented for estimating rapidly the 
aeroelastic effects on the lateral control of swept and unswept 
wings &t subsonic and supersonic speeds. "The charts and 
&pproximate formulas presented herein together with those 
presented in NACA Report 1140 also serve to simplify 
design procedure in many instances because they can be 
used at the preliminary design stage to estimate the amount 
of additional material required to stiffen а wing which 15 
strong enough and because they indicate that the best way 
of distributing this additional material in most cases is to 
locate most of it near the wing tip. 

For the purpose of making specific calculations, the limita- 
tions of the method of this report are that they do not apply 
directly to wings with very low aspect ratio, with very large 
angles of sweep, with zero taper ratio, or with large sources 
of concentrated aerodynamic forces. 
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The charts and, approximate formulas indicate. that the 
control effectiveness of an airplane may be increased by 
varying some of the design parameters such as the ratio of 
torsional to bending stiffness and, if necessary, resorting to 
unconventional lateral-control devices. The charts also in- 
dicate that a wing which is strong enough is most likely to 


be affected by losses in lateral control due to wing flexibility 


if it is to operate at high dynamic pressures, if 16 is thin, if 
it has a large angle of sweepback, if it has an elastic-axis 
location relatively far forward on the chord or & location of 
the center of pressure due to aileron deflection far rearward 
on the chord, or if it is to operate at transonic or high super- 
sonic Mach numbers. 


LANGLEY ÁERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
Lancer Етегр, Va., March 7, 1962. 


APPENDIX 


METHODS OF CALCULATIONS ON WHICH THE CHARTS ARE BASED 


THE AEROELASTIC EQUATIONS 


The assumptions made in the following analysis are the 
same аз those made in reference 2: 

(1) Aerodynamic induction is taken into account by 
applying an overall correction to strip theory. 

(2) Aerodynamic and elastic forces are based upon the 
assumption of small deflections. 

(3) The wing is clamped at the root perpendicular to a 
straight élastic axis, and all deformations are considered to 
be given by the elementary theories of bending and torsion 
about the elastic axis. In addition it is assumed, as in 
reference 1 and elsewhere, that the angle between the 
aileron and the wing is constant along the span of the 
aileron. 

Then, for а wing with an outboard aileron,’ the force per 
unit width on sections perpendicular to the SOT axis 18 


qcC;, | 
l= aq (ob) (А1) 
where 1, is а unit step function of s defined by 
Lals) 
| 
О $ 


5; Sf 


where г; is the spanwise ordinate of the inboard end of the 
aileron, The running torque of this force about the elastic 
axis is 
дс?е Cz, 
aaa (a,;—€a351,) (A2) 
where є is the moment-arm ratio e;/e;. 
The integration of these forces yields the accumulated 
torque and bending moment: 


qc; "eO, 
T= ҮІ ar |" (2 P) (as азд) в 


qc; Cr, ТЫТ с 
“aa ze (72 © lets sds as 


1 


. (A3) 


Men (A4) 


Combining these equations with the equations of elastic 
deformation presented in appendix А. of reference 2 as 


* 


results in two анны. equations of equilibrium: 


d (07 ыл ^ 
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These — аге subject to the following boundary 
conditions: 


“(о cos А--Г sin А)— 06 1.| (А5) 


^c Nle cos A—T sin А)--о481,| (А6) 


= 


Ф(0)=0 (А 7а) 
T(0)—0 (A7b) 
(av 22 i (А7с) 
(ЕТ г „з (Ата) 
ds rns, 
(Ет 2% (Ате) 


The angle of attack due to structural deformations is 


‘related to e and Г by the equation 


a= сов A—T sin A (А8) 


After equations (А5) and (А6) have been solved, the rolling 
moment &bout the wing root may be found from the ex- 
pression 

49 2 Ci 


(A9) 


where the root twisting moment 7' and the root bending 
moment M, are given by equations (АЗ) and (A4) evaluated 
at 8 equal to zero. Then the rolling-moment ratio becomes 


ЖЕТЕЛІ? 
5-1 += (A10) 
үл e D, - В 
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where the functions В, to B, are defined by 


" 1 с ? 0 5 © 
B+eB= | (2) & as (Alle) 
ей f L £O. ast det (A11b) 
с, озб 
1 с 2^ 
B=— f (2) 1, ds* (Allo) 
0 Cr. * 
1 ж. ` 
B= Í | £ 1, ds* ds* (А110) 
0 Јо Cr 5. $ 
and where the parameter d/k is defined by 
d = tan A (A12) 


SOLUTION FOR UNIFORM WINGS 

If the torsional stiffness, the bending stiffness, and the 

chord of the wing have constant values of (GJ),, (Е/),, and 

с,, respectively, опо the wing span, the equations of equilib- 

rium (А5) and (А6) become 

e'cosA— —g*[(e cos А--Г sin A)— ea 1] (А13) 

Г”” вір A=—7 (р сов А-Г sin A) oa 1] (А14) 

where the differentiation denoted by the prime is with re- 
spectto £z jc 
8: 


Differentiating equation (А18) once with respect to £ and 
combining it with equation (A14) yields the single differential 
equation of equilibrium, 


os" доз — q0,— об q*e L/ +9 1, (A15) 
subject to the following boundary conditions: 
a,(1)=0 (А16а) 
o, (0)--0 (A16b) 
ау (0)----47 [0.(0) — e056 1, (0)] (А16с) 


The complete solution of equation (А15) can be readily 
obtained by means of Laplace transforms as 


| 158—256) 50-9 аади 


PUTED к жк ку) 


BON (A17) 


where the functions (2) and f,(£), as well as f,(£) which < 


will be used subsequently, are defined in appendix A of 
reference 2 as 


(В) = Cie- (0, cos vet sin yt) 
(= Се Е eft (о, соз ү? sin yt) 


‚ 9 бее Os сов sin vE) 


where the constants of integration are defined in reference 
2 in terms of the roots of the characteristic equation, —28 
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and B+iy. These three functions are equal to zero when 
Е<0. 

The substitution of this solution into equations (All) 
yields the functions 


Вр, h= ETE -AEH (із) 


Bag", P= D-AIA GH рл) в 
(A18b) 
+ qd LJ — £9 ls 
Bq", k)— KA) JG) — f —£)— ііі +B, (A180) 
В«а*, = 7300) — fs(1 — gH Т А) &)] + 
1 fs) 
ү? а) [fs(1) --) (1 — £)la* uS Lu (A18d) 
Ву= — & (A18e) 
B— &—5 Е? (A18f) 


` The value of 4% at reversal is that value which makes 


С; Cis, in equation (А10) equal to zero for given Ns of 


the parameters К, е, and d. 
For & full-span aileron (£,—1, £,—0), equations (A18) 
become 


(А1да) | 


Виа", у= 91+, 


Bag", =F |де” (Ал) 


* pi LJ) 
Bq", k)= — A (A190) 
Bat, UA) uF) yH BGS ae 
(А194) 
B,——1 (A190) 
H 


and for a “Ир aileron” of very short span (6-0, --0), 
equations (18) become 


Bua" =| A+) kf HAAS *—1} f+ By (A20) 


(A20b) 


BA =| AD 0 Је 


B.(a*,5)— po AEA o 0—16, (209 


Bda ри 1 ла" be (A204) 


(А206) 
(A208) 


Ву= — &, | 
Be= & 
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SOLUTION FOR NONUNIFORM WINGS 


By means of strip theory applied at a finite number of 
points, equation (A1) may be written in matrix notation as 


401, 
UJ ue eel Cos] оо Zo} } (A21) 
and the expression for torque (À2) as 
46102, 
144 * [?1( { as} —e038{ 1.) } (A22) 


The matrix notation used in this analysis is the same as in 
reference 1. 
Equations (A3) and (A4), written in matrix notation, are 


(т}= O M| (EY |е) еә) aam 


(My RE O ИГ (Шағай) (424) 


where the integrating matrices [I’] and [II'] are, respec- 
tively, for single and double integration from tip to root 
апа are given in reference 1. 
matrix station nearest the discontinuity of 1, can be modi- 
fied as in reference 1 in such a manner that premultiplication 
of {1a} by [Г] and [ZJ’,] yields the same area and mo- 
ment about the wing root, respectively, as would be obtained 
by analytic integration. This modification may also be re- 
garded as an attempt to round off the lift distribution near 
the inboard end of the aileron in а physically reasonable 
manner. 

The combination of equations (A23) and (A24) with 
matrix expressions for (о) and {Г} im terms of {T} and 
{M} (see appendix A of ref. 2) yields the equilibrium equation 


[[1]—*LATI а } = — доз { В} (А25) 


which can be solved for { а,}. The aeroelastic matrix [A] is 
defined in reference 1 and the column matrix {8 } is defined 


by 


in Tar in| (c 


«| Чт: | ш [us aso 


After the rolling-moment ratio expressed by equation 
(A10) is set equal to zero, the condition for aileron reversal 
is 


4 (ВВ) Ву Bot eD,4- B,— 0 


where the matrix equivalents of equations (A11) become 


В.+еВ,=2. ГА (| 7% (A278) 
By} eB [IT] H TA (A27b) 
в=—\г]| (2) |а.) (A279 


The value of (1,) at the | 
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Ber =] 1.) (A274) 


Therefore, the condition for aileron reversal expressed in 
matrix notation is 


же (лы 
саб | 2а [ӨЛ UT") Bl (1) (A28) 


Solving equation (A28) for оҙ5 and multiplying the resulting 
equation by (8) yields 


Ее ане 
pal G) Hural аа 


The substitution of this expression for озб{В} into the 
equilibrium equation (A25) yields 


(A29) 


(ора (ААҚ оз) (A30) 
where the aileron-reversal matrix [45] is defined by 
d с с 
(8 ШАП СС) [++ИГ] - 
[Ax] - LA]— E | (2) | Н! (A31) 


ч [чт 


The value of 4% at aileron reversal can be found by the itera- 
tion of this aileron-reversal matrix. 


REPRESENTATION OF RESULTS BY APPROXIMATE FORMULAS 


Approximate formulas, similar to those in reference 2, 
have been used to combine the results of the many compu- 
tations indicated in this analysis; as in the case of the for- 
mulas presented in reference 2, those presented in this sec- 
tion are based on considerations of a semirigid wing. 

The functions В, to B; in equation (А10) have been found 
to be given quite accurately by the following approximate 
expressions: 

— 4*[4*р Қыр 
р с сете 
4% 
И М 4% K+K 
B:=— B, КЕ 4. 1— Kok 


(A328) 


(А821) 
Bie s q"p 
E. q*[q*p (Ag+ Kak 
4 1—Kok 
1 
Ч р 
T 4% to К, 
Bet T 1—К›Ё 
7% 


3— 


(A32c) 


(A32d) 


B;= —B,K; (A3 2e) 


where the factors K, to K, depend on the aileron span, the 
taper ratio, and the spanwise variation of the bending and 
torsion stiffnesses. The factor Ke is independent of the 
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aileron span, and the factor Kg is independent of the stiffness 
variation. 

The following approximate formula may by obtained by 
the substitution of equations (А82) into equation (А10): 


-—Ó 

s s (A33) 
Ci, or 2 
ң zs 


where the value of g at aileron reversal is 


gp (L— Kak) (1 к. a) 
qr= Е (A34) 
(К.К Ead E, + 4+Ки 


With the use of the &pproximate formula presented in 
reference 2, 


Ky 
equation (A34) becomes 
К, (1 -К,$ 4) 
* == (А86) 


1HE --Е.БВ.д«-Ка-Е, ++4-+К = 
The accuracy of equation (А36) compared with the results 
ealculated directly by the method of the preceding sectiori 
is illustrated in figures 6 and 12. 
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Fiaure 12.—Comparison of dynamic йине at aileron reversal 
calculated by the analytic integration method of the appendix with 
those calculated from equation (A36) for swept uniform wings with 
tip ailerons. $8,*—41. д 
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The results of reference 2 indicate that the damping-in-roll 
derivative may be expressed approximately by 


: (A37) 


b C 
iy” -7+ (A38) 


an approximate formula for rolling — T 18 


с ren) 0 


Figure 13 shows the approximate formulas (А88) апа (A39) 
to be in good agreement with more accurately computed 
values. (Actually, a slightly more accurate method of 

estimating C, is given in ref. 2 but the simpler expression | 
given here is comparable in accuracy to the expression for 


к| | | | | | At 
| esgor] | / 
"ANE 


(A39) 


Фо -1.6 -12 -8 
Dynamic - pressure ae yn 


Етаовв 13.— Comparison of rolling power and rolling maneuverability 
calculated by the matrix method of the appendix with those cal- 
culated by equations (А88) and (A39) for uniform wings. 8;*=0.6; 
e— 1.0. 
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С, and its simplicity facilitates its application to the esti- 
mation of the rolling maneuverability.) 


Ап approximate expression for the structural twist due 


to aileron deflection similar to the expressions for the struc- 
tural twists due to geometrical angles of attack given in 
reference 2 has been deduced from the results of the analysis 
in the preceding section: 


о _ 4/Ч9р_ (КаК) Ја КаК Af, 
| 4р 


(А40) 


where f, апа Af, are functions of the spanwise coordinate 
s*, the wing chord and stiffness variations, and the aileron 
span. The accuracy of equation (А40) is indicated in 
figure 14. 
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FIGURE 14.—Comparison of angle-of-attack ratios calculated by the 
matrix method of the appendix with those calculated by equation 
(А40) for uniform wings. 58;*-—0.5; e— 1.0. 
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